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Abstract 

We  suggest  the  M\G\oo  input  process  as  a  viable  model  for  network  traf¬ 
fic  due  to  its  versatility  and  tractability.  To  gauge  its  performance,  we  study 
the  large  buffer  asymptotics  of  a  multiplexer  driven  by  an  M\G\oo  input  pro¬ 
cess.  We  identify  the  process  as  short  or  long-range  dependent  by  means  of 
simple  tests.  The  decay  rate  of  the  tail  probabilities  for  the  buffer  content 
(in  steady-state)  at  the  multiplexer  is  investigated  using  large  deviation  tech¬ 
niques  suggested  by  Duffield  and  O’Connell.  The  appropriate  large  deviations 
scaling  is  found  to  be  related  to  the  forward  recurrence  time  for  the  service 
time  distribution,  and  a  closed-form  expression  is  derived  for  the  correspond¬ 
ing  generalized  limiting  log-moment  generating  function  associated  with  the 
input  process.  Two  very  different  regimes  are  identified.  We  apply  our  re¬ 
sults  to  cases  where  the  service  time  distribution  in  the  M|G|oo  input  model  is 
(i)  Rayleigh  (ii)  Gamma  (iii)  Geometric  (iv)  Weibull  (v)  Log-normal  and  (vi) 

Pareto  -  cases  (v)  and  (vi)  have  recently  been  found  adequate  for  modeling 
packet  traffic  streams  in  certain  networking  applications.  Finally,  we  comment 
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on  the  insufficiency  of  the  short  or  long-range  dependence  in  the  process  in 
clearly  describing  buffer  dynamics. 


1  Introduction 

Several  recent  measurement  studies  have  concluded  that  classical  Poisson-like  traf¬ 
fic  models  are  ill-equipped  to  account  for  time  dependencies  observed  at  multiple 
time  scales  in  a  wide  range  of  networking  applications,  including  Ethernet  LANs 
[12,  16,  26],  VBR  traffic  [13],  WAN  traffic  [25].  As  the  resulting  temporal  corre¬ 
lations  are  expected  to  have  a  significant  impact  on  buffer  engineering  practices, 
this  “failure  of  Poisson  modeling”  has  generated  an  increased  interest  in  a  num¬ 
ber  of  alternative  traffic  models  which  capture  observed  (long-range)  dependencies. 
Proposed  models  include  the  fractional  Brownian  motion  input  model  [20]  and  the 
fractional  Gaussian  noise  input  process  [1];  already  both  have  exposed  clearly  the 
limitations  of  traditional  traffic  models  in  predicting  storage  requirements,  conges¬ 
tion  control  and  other  measures  of  traffic  performance. 

In  this  paper  we  focus  instead  on  the  class  of  M\G\oo  input  processes  as  potential 
traffic  models.  An  M\C\oo  input  process  is  understood  as  the  busy  server  process  of 
a  discrete-time  infinite  server  system  fed  by  a  discrete-time  Poisson  process  of  rate 
A  (customers/slot)  and  with  generic  service  time  cr  distributed  according  to  G.  We 
argue  that  M\G\oo  input  processes  provide  a  viable  alternative  to  existing  traffic 
models;  reasons  which  range  from  flexibility  to  tractability,  are  briefly  presented 
below.  A  lengthier  discussion,  given  later  in  the  body  of  the  paper,  is  based  mostly 
on  the  results  developed  in  the  papers  [23,  24];  we  refer  the  reader  to  these  references 
for  proofs  as  well  as  for  additional  information: 

Firstly,  the  M\G\oo  input  model  has  been  succesfully  investigated  as  a  model  for 
some  wide  area  applications,  e.g.,  Paxson  and  Floyd  [25]  report  a  good  fit  to  TEL¬ 
NET  and  FTP  data  using  a  log-normal  service  time  [25].  However,  the  relevance 
of  the  M\G\oo  input  model  to  network  traffic  modelling  is  perhaps  best  explained 
through  its  connection  to  an  attractive  model  for  aggregate  packet  streams  pro¬ 
posed  by  Likhanov,  Tsybakov  and  Georganas  [18].  They  combine  traffic  generated 
by  several  on-off  sources  with  a  Pareto  distributed  activity  period,  and  show  that 
increasing  the  number  of  sources  yields  a  limiting  behaviour  identical  to  the  M\G\oo 
input  stream  with  a  Pareto  distributed  a.  As  should  be  clear  from  their  analysis, 
the  limiting  result  holds  for  arbitrary  activity  period  distributions,  thereby  provid- 
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ing  a  rationale  for  the  view  that  M\G\oo  input  processes  could  provide  a  natural 
alternative  to  existing  traffic  models,  at  least  for  certain  multiplexed  applications. 
This  limiting  argument  is  similar  to  that  of  using  the  Palm-Khintchin  Theorem  to 
justify  the  Poisson  model  for  interactive  data  traffic. 

Secondly,  the  class  of  M\G\oo  input  processes  has  the  desirable  property  of  being 
stable  under  multiplexing,  i.e.,  the  superposition  of  several  M\G\oo  processes  can 
be  represented  by  an  M\G\oo  input  process. 

Thirdly,  the  M\G\oo  model  is  extremely  versatile  in  that,  dependencies  over  a 
wide  range  of  time  scales  can  be  exhibited  simply  by  controlling  the  tail  behaviour  of 
a  [Prop.  3.1]:  If  P(/i)  denotes  the  autocovariance  of  lag  h  for  the  stationary  version 
of  the  M\G\oo  process,  then 

r(/i)  =  XE[o]e~v\  h  =  0,1,...  (1.1) 

where  —  —  In  P  [a  >  h]  and  a  is  the  forward  recurrence  time  (2.5)  associated  with 
<7.  This  relation  already  indicates  the  tremendous  amount  of  flexibility  in  modeling 
positive  correlation  structures.  The  degree  of  positive  correlation  exhibited  by  an 
M\G\oo  input  process  can  be  further  characterized  by  the  sum  of  the  autocovariances 
(1.1),  or  index  of  dispersion  of  counts  (IDC).  We  show  [Prop.  3.2]  that 

oo  , 

IDC  =  £r(/0  =  -E[a(a  +  l)]  (1-2) 

h= 0  Z 

and  the  process  is  short-range  dependent  (i.e.,  IDC  finite)  if  and  only  if  E  [a2]  is 
finite. 

Temporal  correlations  of  M\G\oo  input  processes  are  known  to  affect  queueing 
performance  [17].  Insights  into  this  phenomenon  can  be  gained  by  analyzing  the 
behaviour  of  a  multiplexer  fed  by  an  M\G\oc  input  process.  For  simplicity,  we 
model  the  multiplexer  as  a  discrete-time  single  server  system  consisting  of  an  infinite 
sized  buffer  and  a  server  with  a  constant  release  rate  c  (cells/slot).  The  number  of 
customers  in  the  input  buffer  at  time  t  is  denoted  by  qt .  Our  performance  index 
is  the  steady-state  buffer  tail  probability  P  [9^  >  b],  as  this  quantity  is  indicative 
of  the  buffer  overflow  probability  in  a  corresponding  finite  buffer  system  with  b 
positions. 

Computing  these  tail  probabilities,  either  analytically  or  numerically,  represents 
a  challenging  problem  in  the  absence  of  any  underlying  Markov  property  for  M\G\oo 
inputs.  Instead,  we  focus  on  the  simpler  task  of  determining  the  tail  behaviour  of 
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the  queue-length  distribution  in  some  asymptotic  sense.  More  precisely,  we  seek 
results  of  the  form 

=  (L3) 

for  some  positive  constant  7*  and  mapping  h  :  JR+  — >  IR+;  these  quantities  are 
characterised  by  A,  G  and  c,  and  should  be  computable  fairly  easily.  In  Section  6 
we  carry  out  the  calculations  for  several  distributions,  thin-tailed  as  well  as  heavy¬ 
tailed  ones.  Drastically  different  behavior  emerge  depending  on  whether  vt  =  0(f) 
or  vt  =  o(t)  [Thms.  5.1  and  5.2]. 

Limits  such  as  (1.3)  provide  a  fair  idea  of  the  tail  of  the  queue-length  distribution, 
and  suggest  approximations  of  the  form 

p  teoo  >  b]  ~  e-/l(6b*  (6  oo).  (1.4) 

Of  course,  the  use  of  the  right  handside  of  (1.4)  to  estimate  P  >  b]  may  be 
fraught  with  difficulties  [5].  Nonetheless,  (1.3)  already  provides  some  qualitative 
insights  into  the  queueing  behavior  at  the  multiplexer,  and  could  in  principle  be 
used  to  produce  guidelines  for  sizing  up  its  buffers. 

Our  focus  here  is  primarily  on  large  deviations  techniques  in  order  to  obtain 
(1.3).  This  approach  has  already  been  adopted  by  a  number  of  authors  [9,  14,  15]. 
Applying  recent  results  by  Duffield  and  O’Connell  [9]  we  can  compute  h(b)  and  7* 
under  reasonably  general  conditions.  Further,  for  a  large  class  of  distributions,  we 
can  select  h(b)  =  ,  and  the  asymptotics  (1.3)— (1.4)  then  take  the  compact  form 

P  [Qoo  >  b]  ~  P  [<7  >  6]7*  ( b  -4  00).  (1.5) 

Hence,  in  many  cases,  including  the  Pareto,  log-normal  and  Weibull  service  times, 
<7oo  and  a  (thus  a)  belong  to  the  same  distributional  class  as  characterized  by  tail 
behavior. 

Sometimes,  in  lieu  of  (1.3),  these  large  deviations  techniques  yield  only  weaker 
asymptotics  in  the  form 

lip  inf  tttt  In  P  [q^  >  6]  >  -7*.  (1.6) 

b—too  n{o) 

This  situation  typically  occurs  when  a  is  heavy-tailed,  in  which  case  large  deviations 
excursions  are  only  one  of  several  causes  for  buffer  exceedances  [15].  Results  such 
as  (1.6)  as  still  useful  in  that  they  provide  bounds  on  decay  rates. 

Comparison  of  (1.5)  with  results  from  [20]  and  [22]  points  already  to  the  complex 
and  subtle  impact  of  (long-range)  dependencies  on  the  tail  probability  P  [ q <*>  >  b ]. 
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Indeed,  in  [20]  the  input  stream  to  the  multiplexer  was  modeled  as  a  fractional  Gaus¬ 
sian  noise  process  exhibiting  long-range  dependence  (in  fact,  self-similarity),  and 
the  buffer  asymptotics  displayed  Weibull-like  characteristics.  On  the  other  hand, 
by  the  results  described  above,  an  M|G|oo  input  process  with  a  Weibull  service  time 
also  yields  Weibull-like  buffer  asymptotics  although  the  input  process  is  now  short- 
range  dependent.  Hence,  the  same  asymptotic  buffer  behavior  can  be  induced  by 
two  vastly  different  input  streams,  one  long-range  dependent  and  the  other  short- 
range  dependent!  To  make  matters  worse,  if  the  pmf  G  were  selected  to  be  Pareto 
instead  of  Weibull,  the  input  process  would  be  long-range  dependent,  in  fact  asymp¬ 
totically  self-similar  [22] ,  but  the  buffer  distribution  would  now  exhibit  Pareto-like 
asymptotics.  To  reiterate  the  main  conclusion  of  [22],  the  value  of  the  Hurst  param¬ 
eter  as  the  sole  indicator  of  long-range  dependence  (via  asymptotic  self-similarity) 
does  not  suffice  for  characterizing  buffer  asymptotics.  Furthermore,  buffer  sizing 
cannot  be  determined  adequately  by  appealing  solely  to  the  short  versus  long-range 
dependence  characterization  of  the  input  model  used,  be  it  of  the  M\G\oo  type  or 
otherwise.  Of  course,  long-range  dependence  (and  its  close  cousin,  self-similarity) 
are  determined  by  second-order  properties  of  the  input  process,  while  aymptotics 
of  the  form  (1.3)  invoke  much  finer  probabilistic  properties.  The  finiteness  of  E  [a2] 
(needed  in  (1.2))  is  obviously  a  poor  marker  for  predicting  the  behavior  of  the 
sequence  {vt,  t  —  1,2,...}  (which  drives  (1.3)).  To  close,  the  diverse  queueing 
behavior,  demonstrated  here  and  tied  to  the  tail  behavior  of  cr,  not  only  confirms 
the  versatility  of  M\G\oo  inputs  as  network  traffic  models,  but  also  points  to  the 
need  for  a  very  careful  and  cautious  approach  in  modeling  network  traffic  when  time 
dependencies  are  either  observed  or  suspected. 

The  paper  is  organized  as  follows:  We  introduce  both  the  class  of  M\GI\oo  inputs 
as  well  as  the  multiplexer  model  in  Section  2  along  with  various  preliminaries.  We 
discuss  the  correlation  structure  of  the  M\GI\oo  input  process  in  Section  3.  Section 
4  develops  general  results  on  buffer  asymptotics  which  are  applied  to  our  specific 
model  in  Section  5.  Finally,  Section  6  illustrates  the  asymptotic  results  for  various 
selections  of  distribution  function  G. 

A  few  words  on  the  notation  used  in  this  paper:  All  rvs  are  defined  on  some 
probability  triple  (0,  J7,  P),  with  E  denoting  the  corresponding  expectation  opera¬ 
tor. 
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2  A  multiplexer  driven  by  M\GI\oo  inputs 


We  model  the  multiplexer  as  a  discrete-time  single  server  queue  with  infinite  buffer 
capacity  which  operates  at  a  constant  rate  and  in  a  first-come  first-served  manner: 
Let  qt  denote  the  number  of  cells  remaining  in  the  buffer  by  the  end  of  slot  [t  —  1,  t), 
and  let  denote  the  number  of  new  cells  which  arrive  at  the  start  of  time  slot 
[t,  t  +  1).  If  the  multiplexer  output  link  can  transmit  c  cells/slot,  then  the  buffer 
content  sequence  {qt,  t  =  0, 1, . . .}  evolves  according  to  Lindley  recursion 

Qo  =  q\  qt+i  =  [qt  +  bt+i  ~  c]+ ,  t  =  0,1,...  (2.1) 

for  some  initial  condition  q. 

Here,  we  account  for  time  dependencies  in  the  cell  input  stream  by  modelling 
the  arrival  process  {bt,  t  =  0, 1, . . .}  as  the  busy  server  process  of  a  discrete-time 
M\G\oo  system.  During  time  slot  [t,  t  +  1),  f3t+ 1  new  customers  arrive  into  the 
system.  Customer  i,  i  =  1, . . .  is  presented  to  its  own  server  and  begins 

service  by  the  start  of  slot  [t  +  l,t  +  2);  its  service  time  has  duration  at+  i,j.  Let 
bt  denote  the  number  of  busy  servers,  or  equivalently  of  customers  still  present  in 
the  system,  at  the  beginning  of  slot  [£,  t  +  1),  with  b  denoting  the  number  of  busy 
servers  initially  present  in  the  system  at  t  —  0. 

The  .CV-valued  rvs  b,  {/3t+i,  t  =  0, 1, . . .}  and  {er^j,  t  =  0, 1, . . . ;  z  =  1,2,...} 
satisfy  the  following  assumptions:  (i)  The  rvs  are  mutually  independent;  (ii)  The 
rvs  {fit+ 1,  t  =  0, 1,...}  are  i.i.d.  Poisson  rvs  with  parameter  A  >  0;  (iii)  The  rvs 
Wt,i,  t  —  1, . . . ;  i  =  1,2, . . .}  are  i.i.d.  with  common  pmf  G  on  {1, 2, . . .}.  We  denote 
by  a  a  generic  W-valued  rv  distributed  according  to  the  pmf  G,  and  throughout  we 
assume  E  [ct]  <  oo. 

No  additional  assumptions  are  made  on  the  rvs  {ao,i,  i  =  1,2, . . .}  which  rep¬ 
resent  the  (residual)  service  durations  of  the  b  customers  present  in  the  system  at 
the  beginning  of  the  slot  [0, 1).  Various  scenarios  can  thus  be  accommodated  within 
this  M\G\oo  model:  If  the  initial  customers  start  their  service  at  time  t  =  0,  then  it 
is  appropriate  to  assume  that  the  rvs  *  =  1,2,.. .}  are  also  i.i.d.  rvs  which  are 
distributed  according  to  the  pmf  G.  On  the  other  hand,  if  we  take  the  viewpoint  that 
the  system  has  been  in  operation  for  some  time,  then  these  rvs  {ao,i,  i  =  1,2, . . .} 
may  be  interpreted  as  the  incomplete  work  (expressed  in  time  slots)  that  the  b 
“initial”  customers  require  from  their  respective  servers  before  their  service  is  com¬ 
pleted.  In  general,  the  statistics  of  the  rvs  {cr0)i,  *  =  1,2,...}  cannot  be  specified  in 
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any  meaningful  way,  except  for  the  situation  that  corresponds  to  the  steady-state 
regime. 

Many  properties  of  M\G\oo  input  processes  derive  from  the  decomposition 

bt  =  b[0) +b[a\  t  =  0,1,...  (2.2) 

where  the  rvs  b[°^  and  b[a^  describe  the  contributions  to  the  number  of  customers  in 
the  system  at  the  beginning  of  slot  [t,t+  1)  from  those  initially  present  (at  t  =  0) 
and  from  the  new  arrivals,  respectively.  It  is  plain  that 

b 

&!0)  =  EiKi>t],  *  =  o,i, (2-3) 

i=l 

and  that  the  rv  b[a^  can  also  be  interpreted  as  the  number  of  busy  servers  in  the 
system  at  the  beginning  of  slot  [f,  t  +  1)  given  that  the  system  was  initially  empty 
(i.e.,  b  =  0). 

In  the  next  proposition  we  state  conditions  for  the  queueing  system  (2.1)  to  admit 
a  steady-state  regime  when  driven  by  the  busy  server  process  {bt,  t  =  0,1,...}. 
Weak  convergence  is  denoted  by 

Proposition  2.1  If  XE  [a]  <  c,  then  there  exists  an  R+ -valued  rv  goo  such  that 
qt  =>t  for  any  choice  of  the  initial  conditions  q,  b  and  {oo i  =  1,2,...}.  The 
system  is  then  said  to  be  stable. 

In  general,  the  busy  server  process  { bt ,  t  =  0,1,...}  is  not  a  (strictly)  stationary 
process,  and  Proposition  2.1  will  not  follow  directly  from  the  well-known  stability  re¬ 
sult  of  Loynes  [19]  for  Lindley  recursions  driven  by  stationary  and  ergodic  sequences. 
The  characterization  of  stability  flows  instead  from  an  extension  of  Loynes’  result  to 
the  case  of  driving  sequences  which  couple  with  their  stationary  and  ergodic  versions 
[2].  To  that  end,  we  show  [21]  that  the  busy  server  process  {bt,  t  =  0, 1, . . .}  indeed 
admits  a  stationary  and  ergodic  version,  thereafter  denoted  by  {tf ,  t  —  0, 1, . . .},  and 
with  which  it  couples  for  any  choice  of  the  initial  conditions  b  and  {cro,*,  z  =  1,2,...}. 
This  stationary  version  {&*,  t  =  0, 1, . . .}  can  represented  through  (2.3)  with 

#(0)  =  £lp»>*],  *  =  0,1,...  (2.4) 

71=1 

where  (i)  the  rvs  {on,  n  =  1,2, . . .}  are  independent  of  the  rv  b  which  is  Poisson 
distributed  with  parameter  AE  [cr],  and  (ii)  the  rvs  {Bn,  n  =  1,2,...}  are  i.i.d.  rvs 


7 


distributed  according  to  the  forward  recurrence  time  a  associated  with  a.  This 
distribution  is  given  by 


gr  =  P[g  =  r]  =  r  =  1)2, . . .  (2.5) 

Throughout,  we  assume  at  a  minimum  the  conditions  of  Proposition  2.1,  and 
write  rin  =  AE  [a]  to  stress  the  fact  that  AE  [<r]  indeed  represents  the  average  input 
rate  into  the  multiplexer.  This  will  come  in  handy  when  comparing  the  effect  on 
the  multiplexer  of  several  traffic  streams  distinguished  by  different  distributions  for 
a  with  a  given  value  for  rin  such  that  rin  <  c. 

3  Correlation  properties 

We  write 

vt  =  -lnP[a  >  t] ,  *  =  1,2, ...  (3.1) 

where  the  forward  recurrence  time  a  associated  with  the  service  time  rv  a,  is  dis¬ 
tributed  according  to  (2.5).  The  following  properties  of  (fo*,  t  =  0, 1, . . .}  are  dis¬ 
cussed  in  [6,  7,  21]. 

Proposition  3.1  The  stationary  and  ergodic  version  {b*,  t  =  0, 1, . . .}  of  the  busy 
server  process  has  the  following  properties: 

1.  For  each  t  =  0, 1, . . .,  the  rv  b*  is  a  Poisson  rv  with  parameter  AE  [a]; 

2.  Its  covariance  structure  is  given  by 

T{h)  =  co v[b*,b}+h] 

=  AE  [(ct  —  h)+] 

=  AE  [a]  e~Vh ,  t,h  =  0,1,2,...  (3.2) 

with  the  convention  vq  =  0. 


Proof.  Fix  h  =  1,2,...:  The  first  expression  for  F(h)  is  well  known  [7].  Next,  note 
that 


m 


AE  [(o  -  h)+] 

OO 

A]CP  [(a~h)+  >r ] 

r=0 
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00 

=  [a  >  h  +  r] 

r=0 

00 

=  A  ^  pk>r] 

r=/i+l 

oo 

=  AEM  Y  p[?  =  r] 

r=h+l 

=  AE[u]P[a>h] 


and  (3.2)  follows  from  (3.1). 


The  strength  of  the  positive  correlation  exhibited  by  the  sequence  {b*,  t  = 
0, 1, . . .}  can  be  formalized  as  follows:  We  say  that  the  sequence  { b *,  t  =  0, 1, ... } 
exhibits  short-range  dependence  if 

OO 

Y  TW  <  oo-  (3-3) 

h=0 

Otherwise,  the  sequence  {&*,  t  =  0, 1, . . .}  is  said  to  be  long-range  dependent  [3,  4]. 

For  M\G\oo  processes  this  dependence  can  be  characterized  through  the  scaling 
{vt,  t  =  1, 2, . . .},  or  alternatively  through  the  finiteness  of  E  [o2] .  First,  from  (3.2) 
we  readily  conclude  that 

Proposition  3.2  Assume  lim4_*oo  ^~t=K.  If  K  <  1  (resp.  >  1),  then  the  station¬ 
ary  sequence  {tf  ,  t  =  0, 1, . . .}  is  short-range  (resp.  long-range)  dependent. 

Proposition  3.3  We  have  the  relation 

OO  \ 

Y  rW  =  AE  M  E  [d]  =  -E  [o(o  +  1)] ,  (3.4) 

h= 0  Z 

so  that  the  stationary  sequence  {6*,  t  =  0, 1, . . .}  is  short-range  (resp.  long-range) 
dependent  if  and  only  if  E  [o2]  is  finite  (resp.  infinite). 


Proof.  From  (3.2)  we  see  that 

OO  OO 

Ym  =  AE[a]^P[a>h] 

h= 0  h= 0 

=  AE  [o]  E  [o] 
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00 

=  AE  [a]  ^  rP  [ct  =  r] 

r= 1 

oo 

=  AE  H  (e  [a])-1  £Vp  [a  >  r] 

r~l 

oo  oo 

=  a|%^p[<7  =  *] 

r=l  t=r 

=  AEP[-4^r) 

t=  1  \r= 1  / 

\  00 

=  ^£*(*+l)Pk  =  *] 

Z  (=1 

and  the  conclusion  (3.4)  is  now  immediate.  ■ 


4  General  Results  on  Buffer  Asymptotics 

Several  authors  [9,  14,  15]  have  derived  asymptotics  such  as  (1.3)  by  means  of  large 
deviations  estimates  associated  with  the  sequence  {St  —  ct,  t  —  0, 1, . . .},  where 

So  =  0;  St  =  b!  +  ...  +  b*t,  t  =  1,2,...  (4.1) 

These  results  have  been  obtained  in  varying  degrees  of  generality,  and  are  summa¬ 
rized  below  as  they  apply  to  the  present  context. 

To  fix  the  terminology,  consider  a  monotone  increasing  1R- valued  sequence  {vt,  t  = 
0, 1, . . .}  such  that  lim  vt  =  oo.  A  sequence  of  R-valued  rvs  {xt,  t  =  0, 1, . . .}  is 

t— >oo 

said  to  satisfy  the  Large  Deviations  Principle  under  scaling  vt  if  there  exists  a  lower- 
semicontinuous  function  1 :  1R  — >  [0,  oo]  such  that  for  every  open  set  G, 

—  inf  I(x)  <  lim  inf —  InP  [xt  6  G]  (4-2) 

xeG  ~  t->oo  vt 

and  for  every  closed  set  F, 

lim  sup  —  InP  [xt  G  F]  <  —  inf  I{x).  (4.3) 

t — >oo  Vt  xe  F 

The  rate  function  I  is  said  to  be  good  if  for  each  r  >  0,  the  level  set  {x  G  1R  : 
I(x)  <  r}  is  a  compact  subset  of  1R. 
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In  many  situations  of  interest,  the  rate  function  can  be  expressed  as  the  Legendre- 
Fenchel  transform  A*  of  another  mapping  A  :  ]R  ->  (— oo,  oo],  namely 


A  *{z)  =  sup{0z  -  A(0)},  z  el R .  (4.4) 


The  reader  is  referred  to  the  monograph  [8]  for  additional  information  on  the  subject 
matter  of  Large  Deviations. 

Consider  two  monotone  increasing  IR+ -valued  sequences  {vt,  t  =  0,1,...}  and 
{at,  t  =  0, 1, . . .}  increasing  at  infinity,  i.e.,  lim^oo  vt  —  limt_,.oo  &t  =  oo,  For  each 
t  =  1, 2, . . .,  we  define 


A t(0)  =  —  InE  exp  (ovt— — —  )  , 
vt  V  at  J  J 


e  e  n. 


(4.5) 


If  we  can  show  the  existence  of  functions  g,  h  :  1R+  — >•  ]R+  such  that  h  is  monotone 

increasing  with  lim  h(b)  =  oo  and  the  limit 
6— >oo 


limM 

6->oo  h{b) 


g{y ),  y  >  o 


(4.6) 


holds,  then  the  following  theorem  obtained  by  Duffield  and  O’Connell  in  [9]  applies. 


Proposition  4.1  Assume  the  arrival  sequence  {bt+ 1,  t  =  0, 1, . . .}  to  be  stationary 
and  ergodic,  and  to  satisfy  the  following  conditions: 

1.  For  each  9  in  JR,  the  limit  A (9)  =  limt_>oo  A t(9)  exists  (possibly  as  an  extended 
real  number); 

2.  The  process  {af1  (St~ct),  t  =  1,2, ...}  satisfies  the  Large  Deviations  Principle 
with  good  rate  function  A*  under  scaling  vt- 

Then,  for  each  y  >  0  we  have 

1K“/®lnP[?”  >  61  -  -9(v)  mf  (4.7) 

If  the  (convex)  rate  function  A*  is  continuous  on  [0,  oo) ,  then  Proposition  4. 1 
immediately  implies  the  lower  bound 

liesfi^)lnP[<,”>',|-_7*  (4-8) 

with  7*  given  by 

7*  =  inf  g(y)A*(y).  (4.9) 
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In  [9],  under  additional  conditions  to  the  ones  of  Proposition  4.1,  a  companion 
upper  bound  to  (4.7)-(4.8)  is  derived.  Proposition  4.2  below  incorporates  these 
conditions  with  a  few  adjustments.  A  different  approach  to  calculating  the  upper 
bound  is  discussed  in  [24],  and  yields  identical  bounds. 


Proposition  4.2  Assume  the  arrival  sequence  {bt+ 1,  t  =  0,1,...}  satisfies  the 
specifications  of  Proposition  4.1  as  well  as  the  following  conditions: 

1.  infx>0  g(x)A*(x)  <  oo; 

2.  For  every  7  >  0,  there  exists  y  =  y( 7)  >  0  such  that 

i. 

j  <  -  “rfp(®)A*(*); 


lim  sup  — — 

&-> 00  h(b) 


In 


V- 


E  ■- 


■7  Vk 


ii. 


iii. 


Then,  we  have 


lim  sup 

6-400 


1 

W) 


lna  1(-)  =  0; 

y 


lim  inf 

t->oo  x>y 


A*(x)vt 

h(xat) 


inf  a(z)A*(x). 

x>y 


lim  sup 

6->oo 


1 

h(b) 


InP  [qoo  >  b]  <  -  inf  g(y)A*(y). 
y>  0 


(4.10) 


When  the  conditions  of  both  Propositions  4.1  and  4.2  are  satisfied,  the  asymp¬ 
totics  distribution  take  the  form  (1.3)  (with  7*  usually  given  by  (4.9)). 

When  the  pmf  G  is  heavy-tailed,  the  conditions  of  Proposition  4.1  hold  whereas 
those  of  Proposition  4.2  do  not,  leaving  us  without  an  upper  bound.  Going  back  to 
the  heuristics  given  in  [15],  we  attribute  this  to  the  fact  that  now  buffer  exceedances 
cannot  be  explained  entirely  by  large  deviations  excursions  in  the  arrival  stream, 
as  there  is  a  need  to  take  into  consideration  the  effect  of  a  single  customer  with  a 
large  workload.  Hence,  any  argument  based  on  large  deviations  techniques  alone 
is  bound  to  fall  short.  However,  we  conjecture  that  (1.3)  still  holds  with  scaling  h 
as  specified  through  (4.6)  but  of  course  with  a  different  value  for  7*.  Work  on  this 
issue  is  in  progress. 


5  Evaluation  of  7*  and  h(b) 

The  key  step  in  applying  the  results  of  the  previous  section  consists  in  finding  two 
monotone  increasing  1R+ -valued  sequences  {vt,  t  =  0, 1, . . .}  and  {at,  t  =  0, 1, . . .} 
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increasing  at  infinity,  such  that  the  limit 


A (9)  =  lim  At(0),  9  e  R  (5.1) 

£-400 

exists  (possibly  as  an  extended  real  number).  It  is  helpful  to  view  at  as  the  scaling 
representative  of  a  law  of  large  numbers  and  vt  as  a  large  deviations  type  scaling. 

In  the  context  of  M\G\oo  processes,  it  is  natural  to  begin  with  the  selection 
at  —  t.  In  the  paper  [23],  we  show  that  the  appropriate  large  deviations  scaling 
{vt,  t  =  1, 2, . . .}  is  given  by  (3.1).  To  state  the  results  more  conveniently,  we  set 


A btt(0)  =  -  InE  [exp(^0St)l  ,  6  e  M 

Vt  t 

(5.2) 

for  each  t  =  1,2,....  Obviously,  if  the  limit 

A b(0)  =  lim  A bt(ff),  9  £  1R 

t— KX) 

(5.3) 

exists,  so  does  (5.1)  with 

A(0)  =  A b(6)  -cO,  0  G  1R 

(5.4) 

and  it  suffices  to  concentrate  on  finding  (5.3).  The  main  facts  along  these  lines  are 
developed  in  the  next  two  theorems;  proofs  are  available  in  [21,  23]. 


Theorem  5.1  Assume  vt  =  0{t)  with  lim^oc  vt/t  =  R  >  0.  Then,  for  each  0  in 
1R,  the  limit  A b{9)  =  lim^oo  A b,t{&)  exists  and  is  given  by 


Ab(0)  = 


AE[ct](^^)e(0)  if  6  <  1 
oo  if  9  >  1 


(5.5) 


where 

X{6)  =  1  + (l-e-Re)  (j^exp  (r(0f?-^)^,  0eIR.  (5.6) 
Moreover,  E(0)  is  finite  for  6  <  1. 


We  say  that  the  sequence  {vt/t,  t  —  1,2,...}  is  monotone  decreasing  in  the  limit 
if  there  exists  a  finite  integer  T  such  that  the  tail  {vt/t,  t  =  T  +  1,T  +  2, . . .}  is 
monotone  decreasing. 


Theorem  5.2  Assume  =  o(t)  with  {vt/t,  t  =  1,2,...}  monotone  decreasing  in 
the  limit.  Assume  further  that  there  exists  a  mapping  T  :  IN  — ¥  M  such  that  (i) 
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r(f)  <  t  for  all  t  -  1,2,...,  (ii)  lim^oo vt^p-  =  oo  and  (iii)  lim^oo  ^ =  0. 

I  (t) 

Then,  for  each  9  in  1R,  the  limit  A b{9)  =  lim^oo  A b,t(9)  exists  and  is  given  by- 


Abie)  = 


AE[ct]0  if  0  <  1 
oo  if  9  >  1. 


(5.7) 


Neither  of  the  above-mentioned  theorems  deals  with  the  case  8—1.  However, 
it  can  be  shown  [21]  that  lim^oo  A^t(l)  does  exist  for  the  examples  discussed  here. 
In  fact,  a  little  thought  indicates  that  the  existence  of  this  limit  suffices  for  our 
purpose,  in  that  its  specific  value  is  not  of  any  consequence  in  evaluating  7*  as  given 
by  (4.9).  Indeed,  in  that  case,  under  the  assumptions  of  either  Theorem  5.1  or  5.2, 
we  see  from  (4.4)  and  (5.4)  that 


A  *(z)  —  sup{0z  —  (A  b(9)  —  c8)} 

8e  R 

=  sup{(c  +  z)8  -  Ab(0)} 

9<  1 

=  sup{(c  +  z)9  -  A6(0)},  z  G  IR  (5.8) 

e<\ 

where  in  the  last  step  we  have  used  the  fact  limtffi  Ah{0)  <  Afc(l).  Moreover,  in  the 
case  vt  =  o{t),  further  computations  are  possible  under  Theorem  5.2:  Prom  (5.8) 
and  the  stability  condition  AE  [cr]  =  rm  <  c,  we  get 


A*(z)  =  sup(z  —  (AE  [ct]  —  c))  9 
e<  1 


=  z  +  c-ri: 


z  >  0, 


whence 


7*  =  inf9(y)(y  +  c-?m). 


y>0 


(5.9) 


(5.10) 


As  becomes  apparent  through  the  examples  discussed  in  [24],  the  term  vt  often 
takes  on  a  form  which  is  computationally  inconvenient  when  checking  the  various 
technical  conditions.  Fortunately,  only  the  asymptotic  behavior  of  vt  matters.  More 
precisely,  consider  another  1R+ -valued  sequence  {wt,  t  —  1,2,..,}  which  is  asymp¬ 
totically  equivalent  to  {vt,  t  —  1,2,...}  in  the  sense  that 


lim  ^  =  1. 

£-><x>  Vt 

In  reference  to  (4.6),  we  now  seek  mappings  h,g  :  1R+  — >■  ]R+  such  that 


6-K50  h{b) 


=  0(y)>  y  >  °. 


(5.11) 


(5.12) 
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(5.13) 


It  is  then  easy  to  check  from  (5.11)  that 

ihn^M=  itaM  hM  =J(I)),  9>0 

b-+oo  h(b)  b— >oo  h(b)  f>-> oo  1C[6/y] 


and  these  mappings  also  satisfy  (4.6)!  A  similar  approach  applies  when  checking 
conditions  (i)-(iii)  in  Theorem  5.2.  Here,  suppose  that  we  have  found  a  mapping 
T  :  W  ->  W  such  that  (i)  T(t)  <  t  for  all  t  —  1,2,...,  (ii)  lim^oo  wt =  oo  and 
(iii)  lim^oo  ^ =  0  -  in  other  words,  conditions  (i)-(iii)  with  respect  to  wt- 
Again,  the  asymptotic  equivalence  (5.11)  implies  that  T  satisfies  conditions  (i)-(iii) 
(with  respect  to  vt)  in  Theorem  5.2. 

In  short,  in  applying  the  results  of  Sections  4  and  5,  we  see  that  the  appropri¬ 
ate  conditions  can  all  be  checked  by  replacing  vt  by  any  other  scaling  wt  which  is 
asymptotically  equivalent  to  it  in  the  sense  of  (5.11),  and  hopefully  more  tractable 
analytically.  We  shall  refer  to  any  such  scaling  as  an  auxiliary  scaling,  and  will  use 
it  to  check  various  conditions. 

If  the  function  u  :  1R.+  — >  R+  is  regularly  varying  and  monotone,  then 


lim 

oo 


u([b/y]) 

u(b) 


=  s(y)> 


y  >  0; 


(5.14) 


in  fact,  it  is  well  known  [11]  that 


g(y )  =  yp,  y>  0;  -oo  <  p  <  oo. 


(5.15) 


This  suggests  a  natural  choice  for  h  as  follows:  Denote  by  w  :  Il+  — >  1R+  the 
piecewise-continuous  inperpolation  of  the  auxiliary  scaling  sequence  {wt,  t  =  0, 1, . . .} 
(with  wo  =  0).  If  w  is  regularly  varying,  then  we  can  select  h(b)  =  w(b )  for  all  b  >  0, 
in  which  case  the  asymptotics  (1.4)  takes  the  form 

P  [qoc  >  b]  ~  P  [S  >  b]1*  (i b  ->  oo).  (5.16) 

A  closer  look  at  (5.16)  indicates  that  for  a  number  of  distributions  the  tail  behaviours 
of  the  service  time  a  and  of  the  queue  length  are  of  the  same  type. 


6  Examples 

We  illustrate  these  asymptotic  results  for  various  choices  of  the  pmf  G;  in  each 
case  we  identify  h(b)  and  compute  7*  in  closed  form  (whenever  possible).  When  a 
simple  expression  is  not  provided  for  7*,  it  should  be  understood  that  its  numerical 
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value  can  obtained  by  solving  an  easy  one-dimensional  optimization  problem.  The 
results  are  presented  in  order  of  increasing  weight  on  the  tail  of  G,  or  equivalently,  in 
order  of  increasing  time  dependence.  We  begin  with  the  Rayleigh  distribution  which 
exhibits  the  least  degree  of  time  dependence  and  proceed  to  discuss  the  Gamma, 
geometric,  Weibull  and  log-normal  cases,  all  of  which  are  short  range  dependent. 
Finally,  we  discuss  the  Pareto  distribution  which  exhibits  long-range  dependence 
(in  fact  asymptotic  self-similarity).  Details  of  the  calculations  are  available  in  [24], 
A  rv  A  is  described  as  having  a  Rayleigh  distribution  with  parameter  a  >  0  if 

P[X<x\  —  l  —  e  2^,  x  >  0.  (6.1) 


The  pmf  G  —  {gr,  r  =  1,2, . . .}  of  the  rv  a  is  said  to  be  an  (integer-valued) 
Rayleigh  distribution  with  parameter  a  >  0,  if  o  =st  [A] ,  whence 


9r  —  C 


(r-l)2 

In'* 


1,2,... 


(6.2) 


Proposition  6.1  If  G  is  a  discrete  Rayleigh  distribution  with  parameter  a  where 
a  >  0,  then 

6^k»  ^2  [Qoo  >  ^]  —  — 7Rayleigh  (6-3) 

for  some  finite  constant  7j*laylcigh . 

A  rv  A  has  a  Gamma  distribution  with  parameters  c  >  0  and  b  >  0  if 

P[A<*]  =  1-^^,  a:>0  (6.4) 

where 

r  oo 

Q(ri,x)=  /  e~Hn~ldt,  r\  >  0,  x  >  0 
Jx 

is  the  incomplete  T-function  and  r(r/)  =  Q(r],  0). 

The  pmf  G  =  {<?r,  r  =  1,2,...}  of  the  rv  o  is  said  to  be  an  (integer-valued) 
Gamma  distribution  with  parameters  c  >  0  and  b  >  0  if  o  —st  [A],  which  yields 

ffr  =  (Q(^,c(r  -  i))  -  Q(6,cr)),  r  =  l,2,...  (6.5) 

Proposition  6.2  IfG  is  a  discrete  Gamma  distribution  with  parameters  c  >  0  and 
b  >  0,  then 

K™  l P  fee  >  b]  =  -TGamma  (6-6) 

where  7Gamma  is  a  finite  constant. 
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The  geometric  pmf  G  =  {gr,  r  =  1,2,...}  with  parameter  0  <  q  <  1,  is  given  by 

gr  =  P  [a  =  r\  =  (1  -  q)qr~1,  r  =  1,2,...  (6.7) 

Proposition  6.3  If  G  is  a  geometric  pmf  of  parameter  q,  with  0  <  q  <  1,  given  by 
(6.7),  then 

lim  -  In  P  [9oo  >  b]  =  -TGeometric  (6-8) 

6->oo  0 

where  Tceometric  JS  a  constant. 


A  rv  X  is  said  to  be  a  Weibull  rv  with  parameters  a  and  /3,  a  >  0  and  0  <  (3  <  1, 


if 


P[X  <x]  =  l-e~ax0,  x  >0 


(6.9) 


The  pmf  G  =  {gr,  r  =  1, 2, . . .}  of  the  rv  a  is  said  to  be  an  (integer-valued)  Weibull 
distribution  with  parameters  a  and  (3  if  o  =st  [X],  in  which  case  we  have 


9r  =  e-<r-V0  -e~ar\  r 


1,2,... 


(6.10) 


Proposition  6.4  If  G  is  a  discrete  Weibull  distribution  with  parameters  a  and  (3, 
a  >  0  and  0  <  f3  <  1,  then 


>  b]  —  “TWeibull 


where 

*  a  (Pic-rin)^0 

TWeibull  —  ^  l  -  (3  ) 


(6.11) 


(6.12) 


A  rv  X  is  said  to  be  a  log-normal  rv  if  X  =st  exp(T)  where  Y  is  a  Gaussian 
rv  with  mean  p  and  variance  62.  The  pmf  G  =  {gr,  r  =  1,2, . . .}  of  the  rv  o  is 
said  to  be  an  (integer-valued)  log-normal  distribution  with  parameters  p  and  6  if 
a  =st  [X].  It  is  easy  to  check  that 


9r 


(6.13) 


where  m  =  e*4,  and  $  is  the  cumulative  distribution  function  of  a  Gaussian  rv  with 
zero  mean  and  unit  variance. 


17 


Proposition  6.5  If  G  is  a  discrete  log-normal  distribution  with  parameters  p  and 
S  as  described  above,  then 


where 


lim  inf 

b— >oo 


1 

(In  b) 2 


InP  [Qcc,  >  b]  >  O'Lognormal 


^Lognormal 


c-rin 
2  S2 


(6.14) 


(6.15) 


In  all  the  cases  considered  so  far,  E  [<r2]  is  finite,  and  the  process  {b* ,  t  =  0, 1, . . .} 
is  therefore  short-range  dependent  by  virtue  of  Proposition  3.3.  However,  in  the  log¬ 
normal  case,  we  have  only  obtained  the  lower  bound.  Going  beyond  the  technical 
conditions,  we  explain  this  departure  from  previous  short-range  dependent  cases 
by  the  following  argument:  Although  the  log-normal  distribution  has  finite  second 
moment,  thereby  insuring  short-range  dependence  (though  barely  as  it  turns  out), 
its  tail  has  become  too  heavy  ( vt  —  o(t ))  to  neglect  the  effect  of  a  single  cutomer 
with  a  large  workload. 

Finally,  the  pmf  G  =  {gr,  r  =  1, 2, . . .}  is  said  to  be  a  Pareto  distribution  with 
parameter  a,  1  <  a  <  2,  if 


lim 

r->  oo 


P  [<7  >  r] 

rp  a 


=  1. 


(6.16) 


For  the  sake  of  concreteness  we  use 


gr  =  P[er  =  r]  =  r  a  -  (r  +  1)  “,  r  =  l,2,...  (6.17) 

so  that  E  [o]  <  oo  while  E  [o2]  =  oo,  and  we  conclude  to  the  long-range  dependence 
of  the  process  {b*,t  =  0, 1, . . .}. 

Proposition  6.6  If  G  is  a  discrete  Pareto  distribution  with  parameter  a,  1  <  a  <  2, 
then 

iiminf  j^lnP  [qoo  >  b]  >  -7pareto  (6-18) 

where 

TPareto  =  («  ~  1)  (c  -  rin)  .  (6.19) 

As  a  final  word  we  recall  our  earlier  comments  at  the  conclusion  of  Section  4; 
when  the  service  time  is  heavy-tailed  as  illustrated  by  the  log-normal  or  Pareto 
cases,  the  conditions  of  Proposition  4.2  are  not  all  satisfied.  The  investigation  of 
the  buffer  asymptotics  for  such  processes  will  require  that  we  look  beyond  the  large 
deviation  techniques  used  thus  far. 


18 


References 


[1]  R.G.  Addie,  M.  Zukerman  and  T.  Neame,  “Fractal  traffic:  Measurements, 
modeling  and  performance  evaluation,”  in  Proceedings  of  Infocom  ’95,  Boston 
(MA),  April  1995,  pp.  985-992. 

[2]  F.  Baccelli  and  P.  Bremaud,  Elements  of  Queueing  Theory:  Palm-Martingale 
Calculus  and  Stochastic  Recurrences,  Applications  of  Mathematics  26, 
Springer- Ver lag,  Berlin  Heidelberger,  1994. 

[3]  J.  Beran,  Statistics  for  Long-Memory  Processes,  Chapman  and  Hall,  New  York 
(NY),  1994. 

[4]  J.  Beran,  R.  Sherman,  M.  S.  Taqqu  and  W.  Willinger  “Long-range  depen¬ 
dence  in  variable  bit-rate  video  traffic,”  IEEE  Transactions  on  Communications 
COM-43  (1995),  pp.  1566-1579. 

[5]  G.L.  Choudhury,  D.M.  Lucantoni  and  W.  Whitt,  “Sqeezing  the  most  out  of 
ATM,”  IEEE  Transactions  on  Communications  COM-44  (1996),  pp.  203-217. 

[6]  D.  R.  Cox,  “Long-Range  Dependence:  A  Review,”  Statistics:  An  Appraisal, 
H.  A.  David  and  H.  T.  David,  Eds.,  The  Iowa  State  University  Press,  Ames 
(IA),  1984,  pp.  55-74. 

[7]  D.  R.  Cox  and  V.  Isham,  Point  Processes,  Chapman  and  Hall,  New  York  (NY), 
1980. 

[8]  A.  Dembo  and  O.  Zeitouni,  Large  Deviation  Techniques  and  Applications,  Jones 
and  Bartlett,  Boston  (MA),  1993. 

[9]  N.  G.  Duffield  and  N.  O’Connell,  “Large  deviations  and  overflow  probabilities 
for  the  general  single  server  queue,  with  applications,”  Mathematical  Proceed¬ 
ings  of  the  Cambridge  Philosophical  Society  118  (1995),  pp.  363-374. 

[10]  A.  Erramilli,  O.  Narayan  and  W.  Willinger,  “Experimental  queuing  analysis 
with  long-range  dependent  packet  traffic,”  IEEE/ACM  Transactions  on  Net¬ 
working  4  (1996),  pp.  209-223. 

[11]  W.  Feller,  An  Introduction  to  Probability  Theory  and  Its  Applications,  Volume 
II,  J.  Wiley  &  Sons,  New  York  (NY),  1966. 


19 


[12]  H.  J.  Fowler  and  W.  E.  Leland,  “Local  area  network  traffic  characteristics,  with 
implications  for  broadband  network  congestion  management,”  IEEE  Journal  on 
Selected  Areas  in  Communications  J  SAC-9  (1991),  pp.  1139-1149. 

[13]  M.  Garrett  and  W.  Willinger,  “Analysis,  modeling  and  generation  of  self¬ 
similar  VBR  video  traffic,”  Proceedings  of  SIGCOMM  ’94,  September  1994, 
pp.  269-280. 

[14]  P.W.  Glynn  and  W.  Whitt,  “Logarithmic  asymptotics  for  steady-state  tail 
probabilities  in  a  single-server  queue,”  Journal  of  Applied  Probability  31 
(1994),  pp. 131-159. 

[15]  G.  Kesidis,  J.  Walrand  and  C.S.  Chang,  “Effective  bandwidths  for  multiclass 
Markov  fluids  and  other  ATM  sources,”  IEEE/ACM  Transactions  on  Network¬ 
ing  1  (1993),  pp.  424-428. 

[16]  W.  Leland,  M.  Taqqu,  W.  Willinger,  and  D.  Wilson,  “On  the  self-similar  nature 
of  ethernet  traffic  (extended  version),”  IEEE/ ACM  Transactions  on  Networking 
2  (1994),  pp.  1-15. 

[17]  M.  Livny,  B.  Melamed  and  A.K.  Tsiolis,  “The  impact  of  autocorrelation  on 
queueing  systems,”  Management  Science  39  (1993),  pp.  322-339. 

[18]  N.  Likhanov,  B.  Tsybakov  and  N.D.  Georganas,  “Analysis  of  an  ATM  buffer 
with  self-similar  (“fractal”)  input  traffic,”  in  Proceedings  of  Infocom  ’95, 
Boston  (MA),  April  1995,  pp.  985-992. 

[19]  R.M.  Loynes,  “The  stability  of  a  queue  with  non-independent  inter-arrival  and 
service  times,”  Proceedings  of  the  Cambridge  Philosophical  Society  58  (1962), 
pp.  497-520. 

[20]  I.  Norros,  “A  storage  model  with  self-similar  input,”  Queueing  Systems  -  The¬ 
ory  &  Applications  16  (1994),  pp.  387-396. 

[21]  M.  Parulekar,  Buffer  Engineering  for  Self-Similar  Traffic,  Ph.D.  Thesis,  Elec¬ 
trical  Engineering  Department,  University  of  Maryland,  College  Park  (MD). 
Expected  December  1996. 

[22]  M.  Parulekar  and  A.M.  Makowski,  “Tail  probabilities  for  a  multiplexer  with 
self-similar  traffic,”  Infocom’96,  April  1996,  San  Francisco  (CA). 


20 


[23]  M.  Parulekar  and  A.M.  Makowski,  “Tail  probabilities  for  M\G\oo  input  pro¬ 
cesses  (I):  Preliminary  asymptotics,”  Queueing  Systems  -  Theory  &  Applica¬ 
tions,  submitted  (1996). 

[24]  M.  Parulekar  and  A.M.  Makowski,  “Tail  probabilities  for  M\G\oo  input  pro¬ 
cesses  (II):  Asymptotics”  in  preparation. 

[25]  V.  Paxson  and  S.  Floyd,  “Wide  area  traffic:  The  failure  of  Poisson  modeling,” 
IEEE/ ACM  Transactions  on  Networking  3  (1993),  pp.  226-244. 

[26]  W.  Willinger,  M.  S.  Taqqu,  W.  E.  Leland  and  D.  V.  Wilson,  “Self-similarity  in 
high-speed  packet  traffic:  Analysis  and  modeling  of  ethernet  traffic  measure¬ 
ments,”  Statistical  Science,  to  appear. 


21 


